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Abstract 



Continuum effects in the weakly bound nuclei close to the drip-line are in- 
vestigated using the analytically soluble Poschl-Teller-Ginocchio potential. 
Pairing correlations are studied within the Hartree-Fock-Bogoliubov method. 
We show that both resonant and non-resonant continuum phase space is ac- 
tive in creating the pairing field. The influence of positive-energy phase space 
is quantified in terms of localizations of states within the nuclear volume. 

PACS numbers: 21.10.Gv,21.10.Pc,21.60.-n,21.60.Cs 
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I. INTRODUCTION 



In the theoretical description of drip-line nuclei, the residual coupling between bound 
states and the continuum is an essential element of the physical situation. Standard de- 
scription of many-fermion systems most often invokes the concept of the Fermi sphere of 
occupied states, and correlations of particles occurring mostly in a narrow zone of the phase 
space around the Fermi surface. It is obvious that whenever the Fermi energy is close to 
zero, which by definition is the case in weakly-bound systems, the zone of correlated states 
must incorporate the phase space of particle continuum. 

Among different types of correlations which are important in nuclei, in drip-line nuclei 
pairing plays a singular role, because the intensity of pairing correlations is a determining 
factor in establishing position of the last bound nucleus. Indeed, in any isotopic chain, the 
lightest unbound nucleus corresponds to an odd-iV system in which the pairing energy A, 
missing in the binding of the odd neutron, equals to the loss of the mean-field binding when 
this neutron is shaken off, see e.g. discussion in Refs. |2|-^|]. Therefore, study of pairing 
correlations in weakly bound systems is since many years one of the leading subjects of 
analyses predicting properties of very neutron rich nuclei. 

After early analyses, see Refs. [^|-[TT|| and references cited and reviewed in Refs. fll"2| , |T3 



the field has now become very active, due to projected developments in radioactive ion beam 
facilities throughout the world. In particular, methods based on using the Skyrme effective 
interaction and on the Hartree-Fock-Bogoliubov (HFB) approximation have been applied to 
both spherical [fl~2| p~6|| and deformed fL7 -fL9[] drip-line nuclei. Spherical drip-line nuclei have 
also been studied by using the Gogny interaction and the HFB method |l2l , |l6| . In addition, 
methods based on the relativistic mean field approach [flC], with pairing correlations included 
within the Hartree-Bogoliubov approximation, have been extensively used to study various 
phenomena in spherical drip-line nuclei [pl| -p8| . Structure of pairs have also been studied by 



solving the two-body problem exactly, including the continuum effects [^9| , and the influence 



of particle resonances on pairing properties have been analyzed in terms of the BCS method 
WO- 



In the framework of the mean-field approximation, pairing correlations are consistently 
included within the variational Hartree-Fock-Bogoliubov (HFB) approach ||32|| . Such an 
approach is appropriate for studying bound nuclei, i.e., systems with negative Fermi energy, 
in which the effects of continuum amount to ensuring that the obtained many-body states 
are stable with respect to virtual excitations of particles and pairs of particles to positive- 
energy phase space. In this study we do not consider scattering problems, for which the wave 
functions are not localized and describe physical situations of particles being scattered off, 
or emitted by nuclei. On the other hand, bound states, however weakly bound, are always 
localized and discrete, and such are the basic features of the approach discussed in the present 
study. A properly executed variational theory, and such is the HFB approach, always yields 
localized bound states. Moreover, when solved in the coordinate representation, the HFB 
method takes fully into account all the mean-field effects of coupling to the continuum 

The HFB equations are relatively easy to solve in the matrix representation, i.e., after 
expanding the quasiparticle wave functions on a suitable basis. When the harmonic oscil- 
lator (HO) basis is used, the method is inappropriate for a description of weakly-bound 
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states, because of an incorrect asymptotic properties of the HO wave functions flzfl . How- 
ever, recently developed methods, which use the so-called transformed harmonic oscillator 
basis [^,0, combine the simplicity of the basis expansion with the corrected asymptotic 
behavior. In this way, large-scale deformed HFB calculations in drip-line nuclei become 
recently possible Another method proposed []35| for a solution of the HFB equations 
relies on using the natural orbitals; although very promising in principle, it has not been 
fully implemented in practice yet. 

Apart from these attempts, most of the HFB calculations in drip-line systems per- 
formed to date use the coordinate representation. Most of these studies have been re- 
stricted to spherical symmetry, for which the coordinate-space problem amounts to solving 
one-dimensional equations. Three-dimensional solutions for deformed nuclei have also been 
obtained by using hybrid methods of solving the HFB equations in the basis of coordinate- 



space Hartree-Fock (HF) wave functions |T7|-[T9|. All such algorithms use wave functions 



approximated on spatial grids of points, and the continuum discretized by suitable large- 
box boundary conditions. Numerical stability and convergence properties of these methods 
have been thoroughly tested, and are considered to be sufficient for the current mean-field 
applications, however, more efficient techniques, using the so-called basis-spline Galerkin 
lattices [B?J] are also being constructed. 

In view of the fact that pairing properties of nuclei near the stability valley are most 
often treated within the BCS approximation, there have been numerous attempts of using 
this simpler method in drip-line nuclei too, see e.g. Refs. |37|- p9| , |30| , pT| . Various tests show 
P ]ll^|40| -f42| , however, that such an approach is often unstable and/or divergent. A mean- 
ingful use of the BCS approximation to the HFB method requires especially careful and 
conscious treatment, which may, in fact, cost more effort than a straightforward application 
of the more involved, but certain HFB method itself. 



Resonance contribution into the HF+BCS equations was recently studied in Refs. 30.31 



This approximate method can take into account the effect for well separated resonances, 
when the resonance energies do not depend strongly on the box size or the cutoff procedure. 



More refined treatments of particle continuum have also recently become available [43,44 



although they have not yet been combined with pairing correlations. On the other hand, 
intuition gained in solving the BCS problem in nuclei tells us that low-energy continuum 
states must more generally contribute to pairing correlations of drip-line nuclei. In the 
present paper we aim at quantifying this influence in the framework of the HFB method. 

The full-scale HFB calculations and/or the analytically soluble models with realistic 
potential shapes are invaluable tools for understanding phenomena associated with the con- 
tinuum in weakly bound systems. In this work, we shall employ the Poschl-Teller-Ginocchio 
(PTG) average potential, which is an extension of the textbook Poschl- Teller potential 
proposed by Ginocchio j|8|. It belongs to a broader class of analytically solvable 



potentials proposed some time ago by Natanzon El. The PTG potential has the main 



features of the nuclear mean field, namely, flat bottom, diffused surface, and asymptotic 
freedom (i.e., it exponentially vanishes at large distances). For this potential, analytical 
solutions are available for wave functions and energies of the single-particle bound states 
and resonances. These advantages make it very useful for applications to nuclear structure 
problems |50|j5T[] . 

In Sec. [TI|, we recall the most important features of PTG potential, in particular those 
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concerning the parametrization of the diffuse surface. Sec. |(T| presents some essentials of 
resonance phenomena and discusses their relation with poles of the S'-matrix. In Sec. [IV|, de- 
tails of the single-particle wave functions and single-particle resonances in the PTG potential 
are recalled and, in particular, properties of wave function corresponding to weakly bound, 
virtual and resonance states are discussed. Sec. [V] contains main results of this work. Here 
results of the HFB calculations with the PTG input are discussed, and the pairing coupling 
to single-particle states in the continuum is analyzed. Finally, Sees. |VJ and |V 11 



summarize 



perspectives of further investigations and the main conclusions we draw from this study, 
respectively. 



II. POSCHL-TELLER-GINOCCHIO POTENTIAL 

The PTG potential depends on four parameters: A which determines its shape and 
diffuseness, v which defines its depth, s (in fm _1 ) which is a scaling factor that can be 
adjusted to obtain a given mean radius of the potential, and a, which allows to account for 
an effective mass of the particle. In the present study, the effective mass is not discussed, 
and therefore we take a=0, for which the radial PTG potential is defined as: 

V PTG (R) = ^(v(R) + c(R)), (2.1) 
2m 

where m is the mass of free particle, v(R) is the central part, 

v(R) = - A. 2 u Lj (v Lj + y 2 ) + ( 2 - 2 ) 

x(l- y 2 ) ( 2 - (7 - AV + 5(1 - AV) 
and c(R) is a (nonstandard) centrifugal barrier, 

c(R) = L(L + 1) (^pj (1 + (A 2 - l)y 2 ). (2.3) 

Dimensionless variable rj is proportional to the radial coordinate R, 

7] = sR, (2.4) 

and y is implicitly defined by the expression: 

1 



^A 2 " 



arctanh(y) + VA 2 — 1 arctanQ/VA 2 — V 



(2.5) 



where < r\ < +oo and, consequently, < y < 1. 

Parameter A can take any positive value. (Potentials with A<1, for which in Eq. ( |2.5|) 
function arctan changes into arctanh, will not be considered in the present study.) In 
principle, parameters A, v, and s can take different values for every value of L and j. 
However, below we use (Lj)-independent values of shape parameter A and scaling factor s. 
On the other hand, as indicated in Eq. Q2.2|) , values of depth parameter v are independently 
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chosen for each L and j, so as to obtain 'reasonable' single-particle spectra, and in particular, 
to simulate the spin-orbit splitting of single-particle levels which is absent in the PTG 
potential. 

Problem of the nonstandard centrifugal barrier (|2.3|) deserves a few words of discussion. 
At the origin (77 — > 0), y decreases as 77, y — > rj, and therefore, 

sHR) _ ^_±I> (2 . 6) 

becomes the standard centrifugal barrier. On the other hand, at large distances (77 — > 00), 
1—y decreases exponentially, i.e., 

c(R) — ,4A2 e -2A 2 (^o) ; (2.7) 

where 

VA 2 - 1 nr. / x 

r/o = — ^2 — arctanvA 2 — 1, (2.8) 

and the nonstandard PTG centrifugal barrier exponentially disappears. Consequently, in 
the PTG potential all partial waves behave asymptotically as the s waves. Therefore, in 
the following analysis we concentrate on the L=0 bound states and resonances for which 
the centrifugal barrier is absent (see discussion in Sees. ^ and 0). In some numerical 
calculations we replace the PTG centrifugal barrier s 2 c(R) (|2.3|) by the physical barrier 
L(L + 1)/-R 2 , i.e., 

Unfortunately, in such the analytical solutions do not exist. 

For small values of A, the PTG potential is wide and diffuse. In Fig. p] (lower panel) 
we present the PTG potentials for L=0 and A=l, 3, and 7, with parameters s and u^j 
(Table |V1I| ) chosen in such a way as to keep the depth and radius of the potential fixed. For 
A = 1 and £=0 (see Fig. |l|), one obtains the Poschl- Teller potential [f|l|f5| which has been 
widely studied, e.g., in the molecular physics. For larger values of A, the PTG potential 
gets steeper and resembles the Woods-Saxon potential. For relatively large values of A and 
relatively small values of the depth parameter, one may find a small barrier at the edge of 
the potential well. Finally, for still larger A, the PTG potential resembles the finite-depth 
square well potential. The PTG potentials presented in Fig. [j] for A=3 and 7 correspond to 
the profiles which are interesting from the point of view of simulating the nuclear mean field 
within a physical range of the diffuseness. 

In the middle panel of Fig. [I] we present the L=4 radial PTG potentials, again for 



A=l, 3, and 7, and parameters s and vij given in Table |VII| . Similarly, the top panel 
(with the corresponding parameters s and u' L j) shows the L=4 potentials for the PTG 
centrifugal barrier s 2 c(R) replaced by the physical centrifugal barrier. One can see that up 
to a little beyond the minimum, the PTG potential reproduces fairly well the potential with 
the physical barrier. However, at larger distances, the PTG centrifugal term disappears too 
quickly. 
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Nature of the energy eigenstates inside the potential well, as well as the nature of scat- 
tering solutions, depend strongly on the shape parameter A. The smaller is the value of A, 
the broader are the single-particle resonances. Below the critical value of A crit = \/2, the 
resonances entirely disappear, i.e., there are no single-particle resonant states in any partial 
wave anymore. 



III. RESONANCES 

In this study we are interested in properties of the single-particle or quasiparticle reso- 
nances in the nuclear average potential. Quantum mechanics deals with numerous types of 
quasi- stationary states. These states may be classified by the singularities of S-matrix or by 
the mechanism of their production. One-particle shape resonances are perhaps the simplest 
long-lived states. The particle is captured via quantum tunneling effect to the strongly at- 
tractive inner region of the potential. In this way, a quasi-stationary state is formed which, 
again through the tunneling effect, may leave this inner region. Similar in nature to the one- 
particle shape resonances are the one-particle virtual states. For these states, the confining 
barrier is absent but the potential exhibits large jump at the potential boundary region 
which, in turn, causes jump in the particle wavelength. As a result, the quasi-stationary 
state is formed, which slowly penetrates from inner to outer region of the potential. The 
one-particle shape resonances and the one-particle virtual states are specifically quantum 
phenomena which have no counterpart in the classical physics. 

Another kind of resonance states, the so-called Feshbach resonances, are formed when 
incoming particle excites many particles in the parent nucleus and is captured in the inter- 
mediate state for which direct decay channels are closed. The decay of this state is then 
proceeding through the series of de-excitations of a parent system to either initial channel 
or to the state of a total system having a lower positive energy | 52f . Detailed microscopic 



theory of these resonances have been worked out in the case, when the coupling of inner exci- 



tation to decay channels is weak [53,54]. This has lead to the formulation of continuum shell 
model (CSM) |54| and to the description, e.g., of giant resonances as quasi-bound A^-particle 
states embedded in the continuum ]^,Q. The CSM was extended recently to the realis- 
tic multi-configurational shell model which describes the coupling between the many body 
wave functions for bound states and the one-particle scattering continuum. This so-called 
Shell Model Embedded in the Continuum (SMEC), allows for a simultaneous description 
not only of the shell model bound states and resonances but also of the radiative capture 
cross-sections [j57],j5£|. 

Another type of many-channel quasi-stationary states may result due to the near- 
threshold singularity ][5^|)(J. This can happen if the overlap of wave function in a given 



channel with other channel wave functions is small, leading to the effective decrease of the 
channel-channel coupling and, hence, to the long lifetime. This resonance mechanism is 
common in the CSM [|54],f)7[]. 

Three-particle, near-threshold long-lived states constitute another class of resonance 
states which is particularly interesting in the context of drip-line physics. In this class, both 
near-threshold virtual states (S'-matrix poles at negative energies on nonphysical sheets) and 
resonance states (S'-matrix poles in the complex plane) can be formed. Long-lived states 
of this class are possible if there exists near-threshold bound, virtual or resonance state in 
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the two-particle subsystem. Multiple transitions of particles between these states of two- 
particle subsystem are leading to the appearance of effective, long-range exchange forces 
in the three-particle system pT|--|63|j. An analogous exchange process is also possible in the 
four-particle systems |)4j] . The above recollection of most common resonance phenomena in 



nuclear physics is by no means exhaustive. 



A. Resonances as poles of the S-matrix 

In the present Section we recall the standard theory of the S-matrix and introduce 
the so-called virtual states, which may appear in the single-particle phase space for small 
energies, and therefore are important for the discussion of pairing correlations in weakly 
bound systems, see Sec. |V} 

Poles of the S-matrix can be located in four different regions of the complex fc-plane, 
corresponding to four regions of the two-sheet complex energy surface p5| (see Fig. 0). 
The position of a pole determines both the behavior of the respective wave function and the 
physical interpretation of the solution. The first region corresponds to the positive imaginary 
A;-axis. The wave functions in this region are normalizable, negative energy solutions of 
the Schrodinger equation and correspond to the bound states of the system. The second 
region is the negative imaginary k-axis. Here solutions of the Schrodinger equation are not 
normalizable (they are exponentially diverging) and, hence, they are not physical. These 
solutions correspond to negative energies on the unphysical sheet of the energy surface and 
they are said to be virtual or antibound (see Fig. ||c). The third region is the sector between 
the positive, real /c-axis and the bisection of the fourth quadrant. Asymptotically, the 
solutions in this region, which correspond to the resonant states with the complex energy, 
are oscillating and exponentially decreasing functions. The imaginary part of the energy 
which is negative in this case, is interpreted as a width of the state. Finally, poles located 
in the remaining region of the complex fc-plane, are also said to be virtual. 

Figure ^a presents different regions of the complex fc-plane where the poles of the S- 
matrix are located, and the corresponding regions on the two sheets of the energy surface 
(Figs. and c). The two sheets are connected along the real positive semiaxis. The arrows 
in Fig. represent the movement of poles, which results from decreasing the depth of the 
potential well. In the general case, the poles corresponding to a bound state and to an 
antibound state move pairwise and cross at a given point (denoted k in Fig. |2|a). For L ^ 0, 
this point is situated at the origin. For L = 0, k can be found lower on the imaginary axis 
and the determination of its position is in general not trivial (see ||65|| ). After the crossing, 
the two poles move pairwise on the lower part of the complex fc-plane. One considers that 
they are associated with the resonance phenomenon, once the pole on the right half-plane has 
crossed the dashed line: 3F£(A;) = —^s(k) (see Fig. ||c)), so that one can interpret its complex 
energy (3?(e) > 0, 3(e) < 0) as the energy and the width of a resonance, respectively. 
Such poles are situated on the unphysical energy sheet but the lower one can influence the 
positive-energy solution on the real positive energy axis. 

Let us now discuss how the poles of the S'-matrix are related to the widths of resonances. 
For that, let us consider the short ranged potential V(R) which tends to zero sufficiently 
fast when R — ► oo. The asymptotic (R — > oo) solution of the Schrodinger equation can be 
written as: 
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V(R -> oo) = A(Jfe) exp(-ikR) + B(Jfe) exp(+iA; J R) 

~ exp(-ikR) + exp(+ifc.R). (3.1) 

We will be interested in the poles of S(k) when A(k) vanishes and we shall consider only 
those poles which are embedded in the fourth quarter (JSt(k) > 0, $s(k) < 0) of the complex 
fc-plane, i.e., those which are associated with the resonance phenomenon. Near the isolated 
i th pole, S(k) can be written as: 

d (In S(k)) = 1_ 

where ki is the complex pole. 

The number of poles in the quarter 3£(/c) > 0, $$(k) < 0, can be found following the 
residue theorem: 

1 / d\nS(k) „ 
N = —f ^j^dk. (3.3) 

Consequently, 



2-KiJ c dk 



The density of states: 



dN . ,d\nS(k) , AX 



-(f)(1)"' 

can be expressed then as follows: 

, , ,dln S(k) . . 

p = (2m)- 1 d( } >- . (3.6) 

Inserting (|3.2|) into ( |3.6| ), it is then easy to see that the level density has a local maximum 
whenever: 

fc = &(Jfei) , z = l,.-.,jV (3.7) 
The value of density at the maximum is: 

and corresponds to the complex energy: 

(hkif . . 

where e t h is the threshold energy. The density peaks have the Lorentzian shape and the 
full-width half-maximum of i th peak is given by: 
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Following this simple example, we assume that the poles of S-matrix on nonphysical 
energy sheets near the real axis correspond to almost all resonance states. Unfortunately, this 
highly plausible assertion remains only a hypothesis because the relation between resonances 
and the S-matrix poles is not determined so rigorously as the correspondence between the 
bound states and the S'-matrix poles on the real axis of the first energy sheet. To affirm 
the correspondence between the observed resonances and the S'-matrix poles on nonphysical 
sheets, certain conditions should be satisfied. First of all, the potential has to be sufficiently 
analytic and has to fall-off sufficiently rapidly at R — > oo, so that the corresponding S'- 
matrix can be safely continued to the unphysical sheets. These conditions are satisfied for 
the PTG potential, though the examples of potentials where the analytic continuation in 
the eigenvalue problem brings around the redundant solutions are known as well BR 57] . In 



general, it is safe to speak about the resonance phenomenon in practical problems if the 
width is not large, i.e., T n Lj/e n Lj < 1 or, in other words, if the distance of the resonance 
pole from the physical region is small. This latter condition, as we shall show in sect. |IV|, is 
never satisfied in the PTG potential for low-lying, near threshold resonances. 



IV. SINGLE-PARTICLE WAVE FUNCTIONS AND RESONANCES IN THE 
POSCHL-TELLER-GINOCCHIO POTENTIAL 

A. Analytical results 



Following Ref. ||48|| , we express solutions of the radial Schrodinger equation with the PTG 



potential ( |2.1| ) as functions of the variable x: 



X 1 /1 A2N„,2> V* > 



(1 + AV 
A 2 )?/ 2 

where y is a function of the radial coordinate R, given by Eqs. ( |2.4|) and fl2.5|) . In the variable 
x, the Schrodinger equation transforms into the Jacobi equation, and its general solution 
can be expressed by means of the hypergeometric function: 

* k L 3 (R) = XL(R)(-^-y (4.2) 

y F ( L + l+(3+u Lj L + l+(3-u Lj l+ 3 l-x\ 



2' 2 



where 



X l(R) = ^[1 + A 2 - (1 - A a x)]i (iy^) t . (4.3) 

For a given complex momentum k, the wave function is specified by two dimensionless 
parameters: 
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(3 



ik 

'7a 2 



and 



"Li 



The bound states occur for parameters (3 defined by: 

P = PnLj, 

where 

3^ 2 



(4.4) 



(4.5) 



(4.6) 



A- 2 PnLj 



2i> + L f-1 



-A 2 K 3 + tt) 



2n + L + - 
2 



(4.7) 



and n > 0. In this case, the hypergeometric function reduces to the Jacobi polynomial, and 



the solution reads 148 



nLj =A/" nij [l + A 2 -(l-A 2 )x]* 

PnLj L+l 



(4.8) 



where Af n Lj is a normalization factor, and the number of bound states is limited by condition 
PnLj > 0, which ensures that the eigenfunctions vanish at the infinity. 
The bound-state energies are given by 



2„2 



€ nLj 



2m 



where the dimensionless eigenenergies are 



£nLj — ^PnLj- 



(4.9) 



(4.10) 



It can be seen from the above equation that the tail of a bound-state wave functions for 

7] -> 00, 

^ nLj (R) oc e -^ A2 ^ o) , (4.11) 

does not explicitly depend on L, i.e., there is no influence of centrifugal barrier on any partial 
wave. 

If A 2 > 2, the resonances will occur for integer n such that: 



n 



>l[A(A 2 -2)-^(v Lj + 1 -)-l] 



(4.12) 



The dimensionless resonance energy is then: 
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£ 



nLj 



(2n + L + l) 2 (A 2 



(4.13) 



and the dimensionless resonance width is: 

( 3 
7 nLi = 4 \2n + L + - 



x 



(2n + L + |) 2 (A 2 -l 



-a 2 :, 1 

2 



2 A 



In the limit: £ n £j — > 0, the resonance width is: 



7nLj 

and, hence, the ratio: •JnLj/ZnLj ~ 
have: 

7nL? 



(0) 
InLj 



4(2n + L + -) 2 , 



(4.14) 



(4.15) 



oo, for all values of parameter A. For small £ n Lj, we 



(o) 

InLj 



£ 



nLj 



£ 



+ 2-4 



nLj 



(0) 
7nij 



(4.16) 



This ratio depends strongly on A. For large (n, L), the quantities £ n Lj and 7 n £j are propor- 
tional to (2n + L) 2 and their ratio is: 



InLj 



£, 



nLj 



(A 2 - If' 2 
' A 2 - 2 ' 



(4.17) 



i.e. ^nLj/£nLj < 1 for A > (10 + 4V/5) 172 . Therefore, in our numerical examples we present 
results for A=3 and 7, which are the values on two sides of the limiting case defined by the 
widths of resonances being equal to resonance energies. 

Solutions in the continuum can be found by analytically continuing the eigenfunctions 
from a discrete negative energy to positive continuous energy (see Eq. (p~2|)). The solutions 
obtained in this way are proportional to hypergeometric functions. Imposing the boundary 
condition that an incoming wave has the momentum k, one can determine the scattering 



function for each angular momentum |48| . Using the asymptotic behavior of the most general 
solution (J4.2I) one obtains the matrix elements of the S*-matrix: 



S(k) 



l) L+1 exp(2/3[A 2 r/ + lnA] 



T[-/3]T[(L + l + /3 + u Lj )/2] 
T[(3}V[(L + l-(3 + i> Lj )/2} 
T[(L + f + /?-i/ L3 -)/2] 
T[{L + \-P-v Lj )/2\ 



(4.18) 



Expression ( |4 . 1 8| ) yields the matrix elements of S'-matrix for the Schrodinger problem 
with the PTG potential without any restrictions, i.e., ( |1.18| ) contains informations about 
all mathematical solutions, both physical and unphysical ones. The poles of the S'-matrix 
in the variable k = isA 2 /3, correspond to the remarkable solutions |65| depending on the 
asymptotic behavior of solutions. 
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B. Wave functions 

In the following we aim at analyzing the influence of weakly bound states and low-energy 
resonances on pairing properties of nuclei near neutron drip lines. Therefore, knowing the 
analytical solutions available for the PTG potential, we chose three sets of parameters such 
that the 3si/2 state is either weakly bound, or virtual, or there is a low-lying resonance in 
the S\/2 channel. These three physical situations can be achieved by fixing parameters A 
and s and shifting the depth parameter v s i. In the specific examples discussed below, we 
use the parameters listed in Table ||, where also the corresponding energies are given. 

Wave functions of the resonant, virtual, and bound 3si/ 2 states are shown in Fig. |3|. 
The resonant wave function is calculated at the energy equal to the real part of the pole 
shown in Table [TI[ Normalizations of the virtual and resonant wave functions are chosen 
so as to match the height of the first maximum of the bound wave function. The three 
wave functions illustrate properties of single-particle phase space in the situation where the 
3si/2 state leaves the realm of bound states, dwells for a short time in a ghost-like zone of 
virtual states, and then reappears as a decent (although very broad) resonance. From the 
point of view of the asymptotic properties, pertinent to the scattering problems, these three 
wave functions are completely different. On the other hand, from the point of view of their 
structure inside the nucleus (see the inset in Fig. they are almost exactly identical. In 
fact, in the scale of the inset, the resonant wave function cannot be distinguished from the 
bound one. 

In Sec. [V]we also study several cases of different positions of the 2d^/2 states. In particu- 
lar, we use two values of the depth parameter v d 3 , which are listed in Table || together with 
the corresponding energies of the 2<i 3 / 2 PTG resonance and virtual states. Whenever the 
spectrum in all partial waves is required, like in the HFB calculations below, we combine 
v d i =4.900 with all the three values of u s i, and u s i =5.034 with all values of u d 3. 

C. Localizations and phase shifts 

Since the pairing potentials depend self-consistently on the pairing densities, they are 
concentrated inside the nucleus, i.e., they are nonzero inside, and go to zero outside the 
nuclear volume. Therefore, the more the continuum wave functions are concentrated inside 
the nucleus, the more they are sensitive to the pairing coupling. In order to quantitatively 
characterize the wave functions from this point of view, we introduce their 'localization' as 
the norm inside the sphere of radius Rl oc , 

L hjj] = [ RLOC \^(r)\ 2 dr, (4.19) 
Jo 

where the radius Rhoc is, for the purpose of the present study, arbitrarily fixed at Rl oc = 
1.5 x R1/2 = 7.578 fm, while R1/2 is the radius where the potential drops to its half value. 
[Note that the volume element 47rr 2 is included in the definition of wave function ip(r).] For 
the bound states, localization is just the probability to find the particle inside the sphere 
i?<i?Loc- For the continuum wave functions which are not normalizable, localization depends 
on the chosen normalization condition. 
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Below we present results for the continuum wave functions normalized in two ways. 
Firstly, we may normalize them to unity inside the box of a large radius. (The value of 
i?box=30fm is chosen for all results obtained in the present study.) Outside the box, the 
wave functions oscillate and have an infinite norm. In this case, the localization is a fraction 
of the probability to find the particle inside the nucleus relative to the probability to find it 
inside the box. Absolute values of such a localization depend, of course, on the radius of the 
box, however, we are only interested in comparing relative localizations of the wave functions 
with different energies. Secondly, we may normalize the continuum wave functions in such 
a way that they all have a common arbitrary amplitude in the asymptotic region. For wave 
functions in which the volume element 4irr 2 is included, as is the case here, their amplitudes 
do not asymptotically depend on r. Again, the absolute values of this localization depend 
on the value of the common amplitude, but the relative values tell us at which energy the 
given wave function is better concentrated inside the nuclear volume. One may chose other 
normalization conditions (to a delta function in energy or in momentum, for example), but 
we did not find any advantage in looking at the corresponding localizations, and we discuss 
here only the two ways described above. 

Fig. [|a shows values of localizations of the Sin PTG continuum states, calculated ana- 
lytically, for the normalization in the box (solid line) and for the normalization to a constant 
amplitude (dashed line). Calculations are performed for the PTG potential for which there 
is a weakly bound 3si/ 2 state (cf. Table |I|), and therefore, the lowest resonance in the 
S1/2 channel appears high in energy. The wide bump in the localization, which is seen near 
40 MeV reflects the known fact that the S'-matrix pole at (39 — 18*) MeV generates continuum 
states which are more localized than the continuum states far from the resonant energies. 
However, values of localizations near the maximum are only a factor of 2-3 larger than those 
far from the maximum, i.e., one cannot a priori expect that the only continuum states which 
can couple to the pairing field are those close to resonances. 

One can see that the localization obtained from normalizing the continuum states in the 
box (solid line) presents numerous wiggles, which appear when the consecutive half-waves 
enter the box. On the other hand, the localization obtained from the amplitude normaliza- 
tion is given by a smooth curve. Both localizations are fairly similar, and therefore, different 
normalization prescriptions do not affect our conclusion about the relative localizations of 
the continuum states. 

There is some difference between these two different normalization prescriptions when 
e— >0. At the origin, (p(r)=A sm(kR). In case of the normalization in the box, we can assume 
for fc— >0 that this expression is valid in the whole box. Then : J(f Loc (j) 2 {r)dr/ J o R 2 (r)<ir ~ 
j^loc r ^^ r j jR r 2^ r _ (i? Loc/ /i?) 3 = const. We have verified this assertion also numerically. In 
the case of normalizing to an amplitude A, for k— >0 one gets: /<f Loc (j) 2 (r)dr ~ \A\ 2 k 2 r 3 /3 — ► 
0. Therefore, the dashed curve turns sharply down when approaching e=0. 

In Fig. |]a we also present localizations of the continuum states calculated numerically 
(circles) in the same box of _Rbox = 30fm, and using the discretization of wave functions on 
a mesh of nodes equally spaced by if=0.25fm. One can see that the numerical results 
perfectly reproduce the analytical calculations for the same normalization. In the numer- 
ical calculations, the box plays merely a role of selecting from the infinite continuous set 
of positive-energy solutions a discrete subset of wave functions which vanish at the box 
boundary. Apart from that, the numerically calculated wave functions are very precise 
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representations of the exact wave functions for some specific discretized values of the energy. 

Fig. |]b shows the phase shifts of the same continuum states calculated analytically 
(solid line) and numerically (circles). (In this case, normalization of the continuum wave 
functions does not play any role.) Again, one can see that the numerical results very precisely 
reproduce the analytical ones in the whole range of studied energies. 

Usually, one identifies the resonance when the phase shift passes tt/2 (or better mr/2, 
as the phase shift is defined modulo 2ir). This definition works well for narrow, well sepa- 
rated resonances. In the case of PTG potential, the resonances are broad and this definition 
is inadequate. It is then better to identify the resonance with the inflection point in the 
derivative d5(e)/de. This also demonstrates real difficulty in identifying resonances in po- 
tentials with diffuse surfaces and proves again the advantage of the soluble models where 
the S-matrix is analytically known and the poles can be analytically studied. 



V. PAIRING IN WEAKLY BOUND SYSTEMS 

A variational mean-field approach to pairing correlations results in the HFB equations 
32fl . In weakly bound systems, these equations should be solved in coordinate space in 
order to properly take into account the closeness of the particle continuum [0,0]. In the 
most general non-local coordinate representation, the HFB equations have the form of the 
following matrix integral eigenequation: 



A/r/^K^') h(ra,r'a') \fc 
jf\h(ra,rW) -h(ra,r'a') J\c 

(E + X \ ( ME,ra)\ 
[ E - X ) \ (j) 2 (E,ra) J 



h(E,r'a') 



(5.1) 



where E is the quasiparticle energy, A is the Fermi energy, and h(ra, r'a') and h(ra, r'a') 
are the mean-field particle-hole (p-h) and particle-particle (p-p) Hamiltonians, respectively. 
Contrary to the HF equations, which define one-component (single-particle) wave functions 
[the eigenstates of h(ra,r'a')}, the HFB method gives two-component (quasiparticle) wave 
functions (the upper and lower components are denoted by <fii{E, ra) and 4>2{E, ra), respec- 
tively). 

In the following, we solve the HFB equations ( |5.1| ) by fixing the p-h Hamiltonian to be 
equal to the sum of the kinetic energy (with constant nucleon mass) and PTG' potential 
( p.9| ). In this way, we study self-consistency only in the pairing channel, while the single- 
particle properties are kept unchanged, and under control. For example, the single-particle 
energies and resonances do not change during the HFB iteration, and are not affected by 
the pairing properties, which would have not been the case had we allowed the usual HFB 
coupling of the p-h and p-p channels. Moreover, within such an approach we only need to 
solve the HFB equations for neutrons, i.e., for the particles which exhibit the weak binding 
under study here. Note that for the physical centrifugal barrier included in all the L>0 
partial waves, the energies of single-particle bound states and resonances are not given by 
analytical expressions (|4.10| ) and fl4.13|) . However, the barrier does not appear in the Si/ 2 



channel, and these energies are still given analytically. 
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Two parameters of the PTG potential have been fixed at values used in the previous 
Sections, namely, A=7 and s=0. 04059, while the depths parameters u^j (Tables |I| and III) 



have been chosen in such a way j5| as to simulate a hypothetical single-particle neutron 
spectrum in drip-line nuclei with iV<82. For the scope of the present study, details of this 
spectrum are insignificant; we only aim at realizing the physical situation where the PTG' 
3si/2 or 2g?3/2 states are near the threshold (close to zero binding energy) and at the same 
time the Fermi energy is negative and small. 

Contrary to the p-h channel, the full self-consistency is required in the p-p channel, with 
the p-p Hamiltonian given by the local pairing potential | 



h{ra, r'a') = U{r)5{r - r% j(T ' , (5.2) 

where 

U(r) = ~V p(r), (5.3) 

and 

P=- E T,MEn,ra)<Pl(E n ,ra). (5.4) 

0<E n <E 

max 



Potential ( |5.3| ) corresponds to the pairing force given by the zero-range interaction, 

V(r 1 -r 2 ) = V 5(r 1 -r 2 ). (5.5) 

Strength parameter Vq has been arbitrarily fixed at Vq=— 175 MeVfm 3 . The pairing phase 
space, given by the cut-off parameter -E max , has been fixed according to the prescription 
formulated in Ref. M. In Eq. (|5.4j) we used the fact that in order to discretize the contin- 



uum HFB states, the HFB equation is solved in a suitable spatial box. The HFB results 
presented below have been obtained with the same box size of i?box=30 fm as those discussed 
in Sec. fVQ . 

As seen from Eqs. ( |5~ll )-( |5~3"l) , the intensity of the pairing coupling [i.e., the off-diagonal 
term in Eq. ( |5.1| )1 is given by the integral of the wave functions with the pairing density 
p(r). This integral can be approximated in the following way: 

J d 3 rJ2<f>*i(En,ra)p(r)<f> 2 (En,r<T) - ~ Po (N n L[4>l]) l/2 , (5.6) 



where £[</>*] is the localization of the upper HFB wave function, defined as in Eq. ( 4.19 ). 
and N n is the norm of the lower HFB wave function: 

N n = f d 3 rY,\MEn,ra)\ 2 . (5.7) 

Equation ( |5.6| ) gives only a very crude approximation, which aims only at showing the 
main trends. It is based on two assumptions: (i) that the pairing potential is constant within 
the radius -Rloc of the sphere for which the localization is defined, and zero otherwise. Zero- 
range pairing force (|5.5f) leads to the volume-type pairing correlations [0, for which the 



pairing densities are spread throughout the nucleus, and can be crudely approximated by a 
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constant value po- Another assumption is: (ii) that the lower and upper HFB wave functions 
are proportional to one another within the radius of -Rloc- We know that this assumption 
holds only in the BCS approximation, while in the HFB approach the lower and upper 
components are different, including different nodal structure 0]. For all energies E n , the 



lower components are localized inside the nucleus, and their norms N n give contributions 



to the particle number (see examples of numerical values presented in Ref. |T2[ and in 
the following subsections). On the other hand, the continuum upper HFB wave functions 
behave asymptotically as plane waves, however, their pairing coupling is dictated by their 
localizations. We are not going to use Eq. ( |5.6| ) in any quantitative way; we only use it 
as a motivation to look at localizations of the upper HFB components as measures of how 
strongly given continuum states contribute to pairing correlations. 



A. The Sij2 continuum 

Figures |] and || show the S1/2 localizations and phase shifts, respectively, of the upper 
HFB quasiparticle wave functions calculated numerically (dots), compared with the analyti- 
cally calculated localizations and phase shifts of the PTG states (lines). We are interested in 
the low-energy si/ 2 continuum, and therefore the Figures show results well below the broad 
(39 — 18*) MeV resonance discussed in Sec. [IV C| . The PTG results (no pairing) are shown 



as functions of the single-particle energy e, while the HFB results are plotted as functions 
of the quasiparticle energy E shifted by the Fermi energy A, i.e., E+X. In this way, at large 
energies the paired and unpaired energy scales coincide. 

Three panels presented in Figs. [5] and |G] correspond to the PTG 3si/2 states being resonant 
(a), virtual (b), and weakly bound (c), cf. Table ||. In order to realize these three different 
physical situations, the bottom of the PTG potential has to be shifted by about 2.75 MeV. 
This illustrates the "width" of the zone where the 3s 1/2 states are virtual. Apparent shift 
in the single-particle energies is much smaller; the 3s±/2 state moves then from the — 74keV 
bound-state energy to the 74keV resonance energy (Table 0). Even if the shift in the 
single-particle states is so small, the corresponding canonical states move by about 1.75 MeV 
(Table [TV]). These latter states appear to be entirely unaffected by the dramatically different 
character of the single-particle 3si/2 states (Fig. |3|). Their wave functions, shown in Fig. [7], 
are almost identical in the three cases (a)-(c). 

Since the canonical states govern the pairing properties of the system, see discussion in 
Ref. P2], their positions explain the changes in the overall pairing intensity (Ajy), and in the 



3si/ 2 occupation t^ an , shown in Table |TV]. It is worth noting that shifting the canonical 3si/ 2 
level from its ~0.8 MeV distance to the Fermi level to a distance of ~2.5 MeV decreases the 
average pairing gap by as much as about 300 keV. 

Apart from the decrease of the pairing intensity described above, there is no other qual- 
itative change in pairing properties when the 3si/ 2 state becomes unbound. The role of 
the weakly bound state is simply taken over by the low-energy continuum. The presence 
and position of the low-energy S1/2 resonance is not essential for the pairing properties. In 
particular, it would have been entirely inappropriate to use this resonance, and its energy 
of 74keV, in any approximate scheme in which the full continuum was replaced by the 
resonances only. 



16 



In all the three cases (a)-(c), in the region of energies between 10 and 20MeV (Fig. [|) 
one can see the "background" localization of the order of 0.2. These energies are far away 
from resonances, and therefore illustrate localizations of a "typical" non-resonant continuum. 
Near the resonances (see Figs. [|a and ||a), localizations are larger, up to 0.3, but do not at 
all dominate over the background value. Therefore, based on estimate ( |5.6| ) one may expect 
that the pairing coupling of the resonant and non-resonant S\/2 continuum is fairly similar. 

The magnitude of this coupling seems to depend primarily on the S\/2 single-particle 
localization in the 2-3 MeV zone above the Fermi surface. (Note that cases (a)-(c) differ 
only by positions of the si/ 2 states; the spectra in other partial waves are kept unchanged.) 
The weakly bound 3si/2 PTG state generates large continuum localization right above the 
threshold, and the opposite is true for the low-lying 3si/ 2 PTG resonance. Therefore, the 
average pairing gap (Ajv) is substantially larger in case (c) than in case (a), and as a 
consequence, the HFB localizations in case (c) differ strongly from the PTG localizations, 
while in case (a) they are almost identical. The same pattern clearly also appears for the 
HFB phase shifts, as compared to the PTG phase shifts, Fig. ||. 

Figure |8] shows norms N n of the lower HFB components, Eq. ( p.7| ). Since the lower HFB 
components are not mutually orthogonal, N n cannot be associated with the occupation 
probabilities. On the other hand, the canonical occupation factors f^ an do play such a role, 
because the canonical states form an orthogonal basis. Comparing values of N n (Fig. |8|) and 
f^ an (Table |TV]), one sees that the canonical 3si/ 2 state collects all the occupation strength 
of the quasiparticle states in the low-energy continuum. Of course, numbers N n scale with 
the overall pairing strength, i.e., they are smaller in case (c) than in case (a), but in every 
case all quasiparticle states below about 5 MeV significantly contribute to f^ an . Moreover, 
choosing only one quasiparticle state in this zone (call it resonance or not), would have 
provided only for about one third of the occupation factor. This illustrates again that the 
Si/2 continuum has to be taken as a real continuum (discretized, if needed), and not through 
any single representative. 

B. The d 3 / 2 continuum 

Since the low-energy si/ 2 resonances are always very broad, the distinction between the 
resonant and non-resonant continuum is not very clear in this channel. On the other hand, 
pairing coupling of very narrow high-j resonances is not very different from that of any other 
bound states, and hence they are not very interesting to look at. Therefore, in addition to 
investigating the Sy 2 resonances, we also study here an intermediate case of low-energy 2<i 3 / 2 
resonances, which can be narrow or broad depending on their positions inside the centrifugal 
barrier. For that we use the PTG' potential which contains the physical centrifugal barrier, 
Eq. fl2.9|) , and chose several different values of the depth parameter v d z. 

First, we study two cases listed in Table [TT], chosen in such a way that the resonances 
are located at two different positions deep inside the centrifugal barrier. (For the chosen 
depths, in the PTG potential the 2c? 3 / 2 states are resonant (d) or virtual (e); their analytical 
energies are given in Table ||.) Although for the physical barrier the analytical results 
are not available, one can estimate the resonance energies to be about e res — (1.5— 0.4i) MeV 
(d) and e res ~(0.9— 0.2i) MeV (c), respectively. In the two cases, the resonances are located 
at about 3.5 MeV and 4.1 MeV below the top of the barrier which is about 5 MeV high. 



17 



One can see that the inclusion of the physical barrier moves the PTG resonance (Table [H]) 
up in energy and decreases its width by a factor of eight, while the PTG virtual state is 
transformed into a true, rather narrow resonance. 

Figure |9] shows the HFB and PTG' localizations of the continuum states. Of course, 
the low-energy PTG' resonances create narrow peaks of localizations at the resonant energies, 
and having resonance widths. Insets in Fig. ^ show that the pairing correlations shift the 
PTG' localizations to slightly higher energies, but otherwise the HFB and PTG' localizations 
are very similar. Beyond the narrow resonances, the HFB and PTG' localizations reach the 
non-resonant background values of about 0.3-0.4, which are almost unaffected by the next, 
very broad <i 3 / 2 resonance at about 40MeV. 

Apart from the change in overall pairing intensity, see Table [V, norms of the lower 
HFB components N n shown in Fig. [10] closely follow the shape of localizations. In order 
to further visualize this property, in Figs. [11] and we show results corresponding to the 
2g?3/2 resonances when they are moved up to the top of the centrifugal barrier. Six panels 
presented in the Figures have been obtained by using the depth parameters z/ d 3 =4.9 (0.1)4.4, 
which corresponds to shifting the bottom of the c/3/2 central potential well from —48.4 to 
— 39.8MeV. As a result, the 2(i 3 / 2 resonances move from about 0.9 MeV to 6MeV (at the 
same time the barrier also slightly increases, from 5 to 6 MeV). Hence, the top panels in 
Figs. |ll| and [12] correspond to a broad resonance located right at the top of the centrifugal 
barrier. 

Localizations of the PTG' c/3/2 continuum states (Fig. [11]) closely follow the pattern 
of broadening and rising resonances. Apart from the lowest two panels, corresponding to 
low-energy resonances, the HFB localizations (dots) do not differ from the PTG' results 
(lines). This is because the larger the distance of the resonance from the Fermi surface 
(Aat~— 0.4keV in all cases), the weaker are, of course, modifications induced by pairing 
correlations. Moreover, the rising resonance leaves behind (at low energies) very low values of 
localizations, which additionally contributes to a decreasing of the overall pairing intensities. 

Interestingly, norms of the lower HFB components N n , presented in Fig. [12], closely follow 
the pattern of localizations. It means that values of localizations of the PTG' single-particle 
states indeed decide about the pairing coupling of the HFB quasiparticle continuum states, 
cf. Eq. (|5.6|) . However, again none of the quasiparticle states (close to resonant energies 
or not) can be used as a single representative of the continuum phase space. Even for the 
resonances located very deeply inside the barrier (the lowest panels in Fig. |l^), norms N n 
for quasiparticle states at resonance energies do not exhaust the occupation numbers v can . 

In Figs. [11] and [12], arrows indicate values of the canonical 2d 3 / 2 energies e can . With 
resonances moving up, these energies increase too, and appear always slightly above the 
resonant energy. These canonical states represent correctly the whole region of the quasi- 
particle continuum within, and above the centrifugal barrier. The 1^3/2 canonical states are 
deeply bound (~— 25 MeV), and the 3<i 3 /2 canonical states are very high (~+25MeV), so 
indeed those corresponding to 2d 3 / 2 (as shown in the Figures) are well separated represen- 
tatives of a wide region of the phase space. Therefore, any approximate scheme, aiming at 
a correct description of pairing correlations near continuum, should vie to find these states 
at right places. Unfortunately, the HF+BCS method which uses resonant continuum states 
only, does not seem to have a chance to attain such a goal. However, similarity of localiza- 
tions and norms N n found in the present study may give hope for approximate solutions of 
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HFB equations. Of course, when the spherical symmetry is conserved, direct solution of the 
HFB problem is easy enough that no approximated methods are necessary. However, this is 
not the case for the deformed systems, where moreover, single-particle resonances split and 
become much more difficult to treat. 



VI. PERSPECTIVES AND OUTLOOK 

The states at positive energies having maximum localization inside the region of potential 
well, preserve a remnant of the eigenvalue structure [^,^]. Whether virtual or resonant, 
these low lying states are likely to play an important role in the low-energy phenomena. 
For example, the neutron direct radiative capture process, which depends on properties 
of nuclear potential both at the surface and in the interior regions, is a unique tool to 
investigate properties of nuclear wave functions. Similar information may also be derived 
from the inverse process, i.e., the Coulomb dissociation process, and both those processes are 
going to be very important in extracting the information about nuclei close to the drip line. 
Raman et al. |7{J have shown how strongly the details of a neutron-nucleus potential affect 



the capture mechanism for s-neutrons. (Extension of this analysis to higher partial waves 
has been done recently by Mengoni et al. J7T| .) This is also important in many astrophysical 
issues, e.g., the inhomogeneous Big Bang model or the neutron poison problem. 

Low-lying resonant states may influence the pairing correlations enhancing the binding 
energy of halo nucleons (neutrons). In this work we have seen that a presence near threshold 
of the L=0 poles of the S'-matrix corresponding to bound, virtual, and resonance states, may 
influence the pairing correlations and increase the localization of scattering wave functions 
in the narrow range of Ae ~2-4MeV above the threshold. Hence, the pairing interaction 
induces a subtle rearrangement effect in the structure of scattering states which modifies 
the strength of pairing field and the binding energy of affected nucleons close to the Fermi 
surface what, in turn induces new rearrangement in the structure of low lying scattering 
states, etc. This effect can be described only in the self-consistent treatment of pairing 
correlations. 

For the HF potentials where the resonance width T tends to a finite value when the 
resonance energy e tends to zero, in the physical situation corresponding to the drip lines, 
i.e., when the nucleon (neutron) separation energy tends to zero (S n — ► 0), and hence 



— A n — A ~0, the standard quasiparticle picture |72| may be questionable. For the quasipar- 
ticle picture to be valid, the appreciable fraction of the single-particle strength should remain 
in the quasiparticle excitations and, moreover, the lifetime of excitations about the Fermi 
surface should be relatively long. Whether the single-particle strength is concentrated in 
the quasiparticle excitations depends on the strength of the residual interaction mixing the 
single-particle strength into more complicated configurations. But even when the residual 
interaction is strong, one can always assure in normal systems that the quasiparticle lifetime 
is long, i.e., the spreading of the quasiparticle strength function is small. In those systems, 
width of a quasiparticle state T(k) ~ (k — /cf) 2 can be made small by letting k approach 
kp |[73|| . This result is independent of details of the interaction and follows from the Pauli 
principle and the phase-space considerations. However, as demonstrated in the example of 
the PTG potential, the states above the Fermi surface of fcp~0 may have a finite width in 
the limit of k — > kp. In spite of the fact that the PTG centrifugal barrier is different from 
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the physical centrifugal barrier, this feature should be seen in the structure of weakly bound 
nuclei near the drip line, at least for systems where L=0 states are present near threshold. 

This particular aspect does not seem to be modified by including the pairing correlations 
within the HFB method. Consequently, the perturbation theory may not be applicable 
in these systems, and various kinds of instabilities caused by the residual correlations are 
expected. These instabilities may change the initial HF(B) vacuum in these nuclei non- 
perturbatively. In other words, one may expect that the spectrum of excitations in the 
unperturbed HF(B) system may not be compatible with the spectrum of excitations in the 
HF(B) system perturbed by the residual correlations. This brings about a new and challeng- 
ing aspect into the experimental studies of shell structure and excitations close to the drip 
line: The atomic nuclei at the drip-line may provide a new kind of quantum open system 
with yet unknown properties of its excitation spectrum, resulting from the strong coupling 
between bound interior states and the environment of scattering states. If discovered, its 
microscopic description may call for new techniques in solving the quantum many-body 
problem. Similar coupling between the localized quasi-stationary states and the scattering 
states has already been recognized to be responsible for the unusual features of the surface 
of nuclear potential |58] . 

One of the most interesting aspects of correlations in drip-line systems is the question 
of multipole instabilities and deformations. A basis to study such phenomena should be 
provided by including the HFB pairing correlations in deformed states. This is a very 



difficult task, and only recently first attempts of such solutions become available [17,18,34 



The problem here is related to solving the HFB equations in a deformed coordinate-space 
representation, in a situation where variational methods are not applicable because the 
HFB equation has a spectrum unbounded from below. However, in our opinion, only such 
an approach may provide a sound basis of quasiparticle states in which other correlations 
(perturbative or not) may further be investigated. 

It is essential to distinguish between five different aspects of the pairing coupling to the 
continuum phase space. First, enough continuum has to be included to cover the zone around 
the Fermi level, which is reasonably larger than the pairing gap. For gaps of about 1 MeV 
the zone which is 5 MeV wide is often used, cf. Ref. ||17|| . Second, the low- lying continuum 
with large localizations has to be taken into account. In cases studied here (Fig. |5]), the 
zone of 5 MeV seems to be enough, however, this aspect strongly depends on the shape and 
depth of the single-particle potential, and the safe limit should probably be at least twice 
larger. Third, enough continuum should be included so as the contributions N n to particle 
number become small, and quasiparticle states entering and leaving this zone do not cause 
significant changes. A limit of iV n ~0. 001-0. 0001 is probably the least one can get away with, 
which already requires going up to 10-15 MeV into the continuum, see Fig. |8|. Fourth, all 
quasiparticle states which significantly contribute to the canonical states (i.e., those which 
have large spectral amplitudes [0) should be included, which requires 10-20 MeV of the 
continuum |L2] . Finally, a coupling between particle-type and hole- type quasiparticle states 
has to be taken into account. Although this coupling is not responsible for the physical 
widths of deep hole states, cf . discussion in Ref. |L2| , it is present in the HFB equations, and 
affects continuum solutions. This aspect requires taking continuum up to energies exceeding 
the depth of the single-particle potential, i.e., to about 40-50 MeV, and such a prescription || 
has been used here. Needless to say that the above discussion concerns the HFB continuum; 
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analyses based on the BCS method, which use the single-particle continuum, have other 
properties, and often diverge with increasing continuum cut-off energy. 



VII. CONCLUSIONS 

Much is known about the analytic properties of one-particle wave functions in the low 
energy continuous spectrum when there is a real, virtual or quasistationary level with energy 
close to zero |68|.|69|]. It was shown for a number of potentials, that the continuum wave 



functions in a wide range of r-values have the r-dependence which is remarkably close to 
that of the wave function for the zero-energy level |7^ , |68|j69[| . Validity of this finding has been 
also discussed in Sec. [IV B| (see Fig. |J). As noticed by Migdal et al. |68| , this approximate 



factorization of the continuum wave functions in the resonance region of one-particle phase 
space should simplify the computation of matrix elements. In the HF+BCS approach of 
Refs. pU|pi| , it is assumed that the analytical features of the S-matrix, which lead to the 



factorization property of the single-particle wave functions in the resonance region, remain 
valid in the presence of pairing correlations. Actually, this weak perturbation hypothesis for 
the S-matrix in the presence of pairing interaction remains not proved. To which extend 
the structure of single-particle resonances and the non-resonant continuum is affected by the 
presence of short-ranged correlations of the pairing type is the main problem which we have 
addressed in this work using the two-component quasiparticle wave-functions of the HFB 
approach and the PTG single-particle potential for which the S'-matrix properties are known 
analytically. The detailed analysis performed in this work, shows that this weak perturbation 
assumption for the analytic structure of the S'-matrix for the one-particle problem may be 
hazardous in many situations. 

The comparison of the localization of HFB upper component of the quasiparticle states 
with the localization of the S1/2 continuum PTG states exposed a complicate interplay 
between resonant and non-resonant HFB continuum which by no means can be approximated 
by the HF+BCS approximation spanned on the skeleton of the S-matrix resonances for the 
one-particle problem. The energy position of canonical states, which govern the pairing 
properties of the system, is not correlated with the presence and position of pole in the 
S-matrix for the one-particle problem. As a consequence, the pairing coupling of resonant 
and non-resonant si/ 2 continuum is quite similar, and its magnitude depends on the single- 
particle localization in the interval of 2-3 MeV above the Fermi surface. 

When the single-particle poles are very close to the real axis in the momentum space 
(§l(k) > 0), like in the case of very narrow high-j resonances deep inside the centrifugal 
barrier, the pairing interaction is too weak to perturb the single-particle pole structure and, 
hence, these resonances are not very different from the bound states. 

In the intermediate case, like <i 3 / 2 -resonances studied in Sec. |VB|, the situation depends 



on the position of the single-particle resonance with respect to the top of the centrifugal 
barrier. In a typical case, however, norms of the lower HFB components closely follow the 
pattern of localizations in the corresponding single-particle continuum, i.e., localizations of 
the single-particle continuum states determine the strength of the pairing coupling of the 
HFB quasiparticle continuum. However, like for the S1/2 continuum, none of the quasiparticle 
states can be used a as a single representative of the continuum phase space. 
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TABLES 



TABLE I. Parameters used to draw the potentials presented in Fig. H. Parameters u 9 and 

v 1 9 are used with the PTG and PTG' centrifugal barriers, respectively. 

9 2 



A 


slfm- 1 } 


b 2 


V 9\ 


9% 


1 


0.2192 


6.719 


6.188 


7.199 


3 


0.09259 


5.165 


5.470 


5.176 


7 


0.04059 


5.034 


5.372 


4.693 



TABLE II. Depth parameters of the PTG potential, for which the 2>S\/2 state is resonant, 
virtual, or bound, and for which the 2d 3 / 2 state is resonant or virtual. The other two parameters 
are A=7 and s=0.04059. 



case 






energy 


(a) 


resonant 3si/ 2 state 


4.882 


e res =(74-4225i) keV 


(b) 


virtual 3s]/ 2 state 


4.972 


e v irt=-84keV 


(c) 


bound 2>Si/2 state 


5.034 


ebound=-74keV 


(d) 


resonant 2d 3 / 2 state 


4.850 


e res =(649-2610i) keV 


(c) 


virtual 2<i 3/ / 2 state 


4.900 


e virt =(-247-1776i) keV 



TABLE III. Depth parameters of the PTG' potential (2.9) defining the neutron spectrum used 
in the HFB calculations. Parameters not given here are listed in Table 0. 



Lj 




Lj 


VLj 


Lj 




Lj 




Sl/2 


Table 


Pl/2 


4.640 


P3/2 


4.880 


d-3/2 


Table 


db/2 


5.180 


fb/2 


4.300 


h/2 


4.720 


97/2 


4.420 


99/2 


4.800 


h-9/2 


4.493 


hn/2 


4.992 


hl/2 


4.493 


*13/2 


4.493 


J13/2 


4.493 
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TABLE IV. Properties of the HFB solutions obtained for the PTG' potentials which give the 
resonant (a), virtual (b), and weakly bound (c) 3si/ 2 single-particle states, or the 2d 3 / 2 low-energy 
narrow resonances at two different positions. The neutron Fermi energies A at and average pairing 
gaps (Ajv) H are given together with the canonical energies e can [12] and occupation factors v 2 an . 
All energies are in keV. 



case 


Xn 


(Aiv) 


Lj 


Scan 


v 2 

u can 


(a) 


-314 


1146 


3^1/2 


2148 


0.0321 


(b) 


-384 


1292 


1043 


0.1002 


(c) 


-488 


1421 




390 


0.2083 


(d) 


-436 


1217 


2^3/2 


1765 


0.0548 


(c) 


-488 


1421 


1107 


0.1149 
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FIGURES 




FIG. 1. Radial PTG potentials for three values of the parameter A (see Table VII for the 
complete set of parameters). The bottom panel shows the potential for L=0, while the middle 
and top panels correspond to L=4 and the PTG fl2,l|) and PTG (2.9) potentials, respectively (see 
text). 
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FIG. 2. Schematic representation of different domains of S'-matrix poles in the complex 
momentum and energy planes. 
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FIG. 3. Radial PTG wave functions of the 3s]/2 state calculated analytically for parameters 
listed in Table ||. The dotted, dashed, and solid lines correspond to resonant, virtual, and bound 
3s 1/2 states, respectively. In the inset the same wave functions are shown in the region of attractive 
potential. Volume element Airr 2 is included in the definition of wave functions. 
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FIG. 4. Localization [Eq. ( [4.19| )] (a) and phase shift (b) of the S1/2 continuum PTG states. 
Circles denote results of numerical calculations performed in the box of 30 fm. Solid and dashed 
lines show the localizations calculated analytically for scattering wave functions normalized in the 
box of 30 fm, and to a constant amplitude of the outgoing wave, respectively. Arrows indicate 
position of the lowest, very broad, resonance. Parameters of the PTG potential are given in 
Table || (bound 3si/ 2 state). 
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FIG. 5. Localization [Eq. ( |4.19 )] of the HFB upper components (dots) of the quasiparticle 
states, compared with the localization of the s^m continuum PTG states (lines). 



32 



T3 
(0 



(0 

(D 
</) 
(0 



3 

2 

1 - 

- 

-1 - 

-2 - 



I I I I | I I I I | I I I I | I I I I 

e res =(74-4225i) keV 



(a) 



3 - 
2 - 




(b) 



e virt=" 84 keV 




■2 - 



5 10 15 

Energy (MeV) 



20 



FIG. 6. Phase shifts of the HFB upper components (dots) of the quasiparticle states, compared 
with the phase shifts of the Si/ 2 continuum PTG states (lines). 
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FIG. 7. Canonical wave functions of the 3sx/2 state, corresponding to the three cases of the 
HFB calculations listed in Table IV. The dotted, dashed, and solid lines correspond to resonant, 
virtual, and bound PTG 3si/2 states, respectively. In the inset the same wave functions are shown 
in the logarithmic scale. Volume element Anr 2 is included in the definition of wave functions. 
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FIG. 8. Norms N n of the S]/2 lower quasiparticle wave functions 4>2(E n ,ra), corresponding 
to the three cases of the HFB calculations listed in Table IV. Inset shows the same values in the 
logarithmic scale. 
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FIG. 9. Localization of the HFB upper components (dots) of the quasiparticle states compared 
with the localization of the d 3 / 2 continuum PTG' states (lines). Upper (a) and lower (b) panels 
differ by the position of the 2d 3 / 2 resonances with respect to the top of the centrifugal barrier, see 
text. 
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FIG. 10. Norms N n of the d%/2 lower quasiparticle wave functions (p2(E n ,ra), corresponding 
to the two cases of the HFB calculations listed in Table IV. Inset shows the same values in the 
logarithmic scale. 
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Energy (MeV) 

FIG. 11. Localizations of the HFB upper components of the d 3 / 2 quasiparticle states (dots) 
compared with those of the d%j2 continuum PTG' states (lines). The six panels differ by the 
positions of the lowest d^/2 resonance with respect to the top of the centrifugal barrier, see text. 
Arrows indicate values of the canonical energies £ can of the 2,d$/2 canonical states. 
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FIG. 12. Norms N n of the HFB lower components (j>2(E n , ra) of the d 3 / 2 quasiparticle states. 
In each panel (similarly as in Fig. [ll]) , arrows indicate values of the canonical energies £can of the 
2^3/2 canonical states, while the corresponding occupation factors v can are given explicitly. 
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